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Cosmic ray antiprotons provide an important probe for the study of cosmic ray propagation in the interstellar 
space and to investigate the existence of Galactic dark matter. 

The ARGO-YBJ experiment is observing the Moon shadow with high statistical significance at an energy 
threshold of a few hundred GeV. Using all the data collected until November 2009, we set two upper limits on 
the p/p flux ratio: 5% at an energy of 1.4 TeV and 6% at 5 TeV with a confidence level of 90%. In the few-TeV 
range the ARGO-YBJ results are the lowest available, useful to constrain models for antiproton production in 
antimatter domains. 



1. INTRODUCTION 

Very High Energy Cosmic Ray (VHE CR) an- 
tiprotons are an essential diagnostic tool to ap- 
proach the solution of several big topics of cos- 
mology, astrophysics and particle physics, besides 
for studying fundamental properties of the CR 
sources and propagation medium. The enigma 
of the matter/anti- matter asymmetry in the lo- 
cal Universe, that of the existence of antimat- 
ter regions, the signatures of physics beyond the 
standard model of particles and fields, as well as 
the determination of the essential features of CR 
propagation in the insterstellar medium, these are 
only a few research topics which would greatly 
benefit from the detection of VHE antiprotons 

BEEEEI. 

Cosmic rays are hampered by the Moon, there- 
fore a deficit of CRs in its direction is expected 
(the so-called 'Moon shadow'). Moreover, the 
Earth-Moon system acts as a magnetic spectrom- 
eter. In fact, due to the GeoMagnetic Field 
(GMF) the Moon shadow shifts westward by an 
amount depending on the primary CR energy. 
The paths of primary antiprotons are therefore 
deflected in the opposite sense in their way to the 
Earth. This effect allows, in principle, the search 
for antiparticles in the opposite direction of the 



observed Moon shadow. 

If the energy is low enough and the angular res- 
olution good we can distinguish the two shadows, 
one shifted towards West due to the protons and 
the other shifted towards East due to the antipro- 
tons. At high energies (> 10 TeV) the magnetic 
deflection is too small compared to the angular 
resolution and the two shadows cannot be disen- 
tangled. At low energies (<500 GeV) the shadows 
are much deflected but washed out by the poor 
angular resolution, thus the sensitivity is limited. 
Therefore, there is an optimal energy window for 
the measurement of the antiproton abundance. 

The ARGO-YBJ experiment is especially re- 
commended for the measurement of the CR anti- 
matter content via the observation of the Galac- 
tic CR shadowing effect due to: (1) good angular 
resolution and pointing accuracy and their long 
term stability; (2) low energy threshold; (3) real 
sensitivity to the GMF. Indeed, the low energy 
threshold of the detector allows the observation 
of the shadowing effect with a sensitivity of about 
9 standard deviations (s.d.) per month at ~TeV 
energy. 

In this paper we report the measurement of the 
p/p ratio in the TeV energy region with all the 
data collected during the period from July 2006 
to November 2009. 
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2. The ARGO-YBJ experiment 

The ARGO-YBJ experiment, located at the 
YangBaJing Cosmic Ray Laboratory (Tibet, P.R. 
China, 4300 m a.s.L, 606 g/cm 2 ), is currently the 
only air shower array exploiting the full coverage 
approach at high altitude, with the aim of study- 
ing the cosmic radiation at an energy threshold 
of a few hundred GeV. 

The detector is composed of a central carpet 
large ~74 x 78 m 2 , made of a single layer of Resis- 
tive Plate Chambers (RPCs) with ^93% of active 
area, enclosed by a guard ring partially (~20%) 
instrumented up to ~100xll0 m 2 . The appara- 
tus has modular structure, the basic data acquisi- 
tion element being a cluster (5.7x7.6 m 2 ), made 
of 12 RPCs (2.8x1.25 m 2 each). Each chamber is 
read by 80 external strips of 6.75x61.8 cm 2 (the 
spatial pixels), logically organized in 10 indepen- 
dent pads of 55.6x61.8 cm 2 which represent the 
time pixels of the detector. The read-out of 18360 
pads and 146880 strips are the experimental out- 
put of the detector [6]. The RPCs are operated in 
streamer mode by using a gas mixture (Ar 15%, 
Isobutane 10%, TetraFluoroEthane 75%) for high 
altitude operation. The high voltage set at 7.2 
kV ensures an overall efficiency of about 96% [7]. 
The central carpet contains 130 clusters (here- 
after, ARGO-130) and the full detector is com- 
posed of 153 clusters for a total active surface 
of ^6700 m 2 . A simple, yet powerful, electronic 
logic has been implemented to build an inclusive 
trigger. This logic is based on a time correlation 
between the pad signals depending on their rela- 
tive distance. In this way, all the shower events 
giving a number of fired pads N pa d > N fr .; g in the 
central carpet in a time window of 420 ns gener- 
ate the trigger. This can work with high efficiency 
down to N tr i g = 20, keeping the rate of random 
coincidences negligible. 

The whole system, in smooth data taking since 
July 2006 firstly with ARGO-130, is in stable data 
taking with the full apparatus of 153 clusters since 
November 2007 with the trigger condition N tr i g 
= 20 and a duty cycle >85%. The trigger rate 
is ^3.6 kHz with a dead time of 4%. The main 
results after about 3 years of stable operation are 
summarized in [5]. 



The reconstruction of the shower parameters is 
split into the following steps. First the shower 
core position is derived with the Maximum Like- 
lihood method from the lateral density distribu- 
tion of the secondary particles [9]. In the second 
step, given the core position, the shower axis is re- 
constructed by means of an iterative un- weighted 
planar fit able to reject the time values belong- 
ing to the non-gaussian tails of the arrival time 
distributions. Finally, a conical correction with a 
slope fixed to a = 0.03 rad is applied to the sur- 
viving hits in order to improve the angular reso- 
lution |10j . Unlike the information on the plane 
surface, the conical correction is obtained via a 
weighted fit which lowers the contributions from 
delayed secondary particles, not belonging to the 
shower front. In detail, we firstly fit a conical sur- 
face to the shower image, by minimizing the sum 
of the squares of the time residual distribution. 
At this stage, all the particles hitting the detec- 
tor have the same weight w»=l. After computing 
the RMS of the time residuals with respect to 
such a conical surface, we set K = 2.5 RMS as 
'scale parameter' and perform the minimization 
of the square of the time residuals weighted sum, 
where w,=l if the particle is onward the shower 
front, Wi=/((i^ p — t{ lt )/K) otherwise. f(x) is 
a common Tukey biweight function. The fit pro- 
cedure is iterated, refreshing every time the scale 
parameter, until the last reconstructed direction 
differs from the previous one less than 0.1°. 

The analysis reported in this paper refers to 
events collected after the following selections: (1) 
more than 20 strips on the ARGO-130 carpet; 
(2) zenith angle of the shower arrival direction 
less than 50°; (3) reconstructed core position in- 
side a 150x150 m 2 area centered on the detec- 
tor. Moreover a quality cut to remove misrecon- 
structed events is applied. According to the simu- 
lation, the median energy of the selected protons 
is E 50 ssl.8 TeV (mode energy «0.7 TeV). 

3. Monte Carlo simulation 

A detailed Monte Carlo (MC) simulation of the 
CR propagation in the Earth-Moon system [11] 
has been developed in order to estimate the ex- 
pected antiproton flux in the side opposite to the 
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CR Moon shadow. The simulation is based on 
the real data acquisition time. The Moon posi- 
tion has been computed at fixed times, starting 
from July 2006 up to November 2009. Such in- 
stants are distant 30 seconds each other. For each 
time, after checking that the data acquisition was 
effectively running and the Moon was in the field 
of view, primaries were generated with arrival di- 
rections sampled within the Moon disc basing on 
the effective exposure time. 

After accounting for the arrival direction cor- 
rection due to the magnetic bending effect, the 
air shower development in the atmosphere has 
been generated by the CORSIKA v. 6.500 code 
with QGS JET /GHEISHA models [H]. CR spec- 
tra have been simulated in the energy range from 
10 GeV to 1 PeV following the results given in 
|13) . About 10 8 showers have been sampled in 
the zenith angle interval 0-60 degrees. The sec- 
ondary particles have been propagated down to 
cut-off energies of 1 MeV (electromagnetic com- 
ponent) and 100 MeV (muons and hadrons). The 
experimental conditions (trigger logic, time res- 
olution, electronic noises, etc.) have been re- 
produced by a GEANT4-based code [T5]. The 
core positions have been randomly sampled in 
an energy-dependent area large up to 2-10 3 x 2- 
10 3 m 2 , centered on the detector. The simulated 
events have been generated in the same format 
used for the experimental data and analyzed with 
the same reconstruction code. 

4. Data analysis 

For the analysis of the shadowing effect, the 
signal is collected within a 10°xl0° sky region 
centered on the Moon position. We used celes- 
tial coordinates (right ascension and declination, 
hereafter R.A. and DEC.) to produce the event 
and background sky maps, with 0.1°x0.1° bin 
size. Finally, after a smoothing procedure, the 
significance map, used to estimate the statistical 
significance of the observation, is built. 

As can be noticed from Fig. [TJ the Moon 
shadow turns out to be a lack in the smooth 
CR signal observed by ARGO-YBJ, even without 
subtracting the background contribution. The 
background events are not uniformly distributed 
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Figure 1. Deficit of CRs around the Moon posi- 
tion projected along the R.A. direction. Showers 
with N>60 recorded from July 2006 until Novem- 
ber 2009 are shown. 



around the Moon, because of the non-uniform ex- 
posure of the map bins to CR radiation. 

The background has been estimated with two 
different approaches, the time swapping and the 
equi-zenith angle methods, as described in [15 , 
in order to investigate possible systematic un- 
certainties in the background calculation. The 
significance map is obtained from the event and 
background maps after applying a smoothing pro- 
cedure to take into account the angular resolution 
of the detector [15] . A detailed study of the two 
background calculation methods has shown that 
they give results consistent with each other within 
1 s.d.. 

A significance map of the Moon region is shown 
in Fig. [2] It contains all the events belong- 
ing to the lowest multiplicity bin investigated 
(20<N<40), collected by ARGO-YBJ during the 
period July 2006 - November 2009 (about 3200 
hours on-source in total). The significance of the 
maximum deficit is about 22 s.d.. The observed 
westward displacement of the Moon shadow by 
about 1.5° allows to appreciate the sensitivity of 
the ARGO-YBJ experiment to the GMF. In such 
a map a potential antiproton signal is expected 
eastward within 1.5° from the actual Moon po- 
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Figure 2. Moon shadow significance map. It col- 
lects all the events detected by ARGO-YBJ from 
July 2006 until November 2009. The event mul- 
tiplicity is 20<N<40 and zenith angle 9 < 50°. 
The color scale gives the statistical significance. 



Figure 3. Upper (lower) panel: deficit counts 
observed around the Moon projected along the 
North-South (East-West) axis for 40<N<100 
(black squares) compared to the MC simulation 
expectations. 



sition (i.e., within 3° from the observed Moon 
shadow). The median energy of the selected 
events is E 50 «750 GeV (mode energy ss 550 
GeV) for proton-induced showers. The corre- 
sponding angular resolution is ~1.9°. However, 
the large displacement of the shadow is only one 
ingredient of this analysis, the other being the 
angular resolution which is not small enough in 
this multiplicity range. Indeed, as can be seen 
from Fig. [5J the matter shadow extends to the 
antimatter side with a significance of about 10 
s.d., thus limiting the sensitivity to the antipro- 
ton abundance measurement. 

We stress that this is the first time that an 
EAS array is able to detect the Moon shadow 
so shifted, observing the signal due to sub-TeV 
primary CRs. 

5. Results and discussion 

The optimal energy windows for the mea- 
surement of the antiproton abundance in CRs 
are given by the following multiplicity ranges: 
40<N<100 and N>100. In the former bin the 



statistical significance of Moon shadow observa- 
tion is 34 s.d., the measured angular resolution 
is ~1°, the proton median energy is 1.4 TeV and 
the number of deficit events about 183000. In the 
latter multiplicity bin the significance is 55 s.d., 
the measured angular resolution ^0.6°, the pro- 
ton median energy is 5 TeV and the number of 
deficit events about 46500. 

The chance of unfolding all CR spectral contri- 
butions relies on MC simulations, as well as the 
search for antiprotons demands to properly re- 
produce the Moon shadow signal. The projection 
along the North-South (East- West) direction of 
the deficit counts around the Moon is shown in 
the upper (lower) panel of Fig. 03 for 40<N<100. 
The vertical axis reports the events contained 
in the angular slice parallel to the North-South 
(East- West) axis and centered to the observed 
Moon position. The width of this band is ±3°. 
The data are in good agreement with the MC sim- 
ulation and the observed shadow is shifted west- 
ward of about 1° (lower panel). 

In order to evaluate the p/p ratio, the projec- 
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tions of the Moon shadow along the R.A. direc- 
tion have been considered. The GMF shifts west- 
ward the dip of the signal from positively charged 
primaries. Searching antiprotons means looking 
for excesses in the eastern part of the R.A. pro- 
jection, i.e. trying to fit the Moon shape expected 
from combining CR and antimatter to the shape 
obtained from the experimental data. It is worth 
noticing that for the antiprotons we assume the 
same energy spectrum of the protons. Of course, 
whichever matter-antiprotons combination is ob- 
tained, the total amount of triggered events must 
not be changed, so that the fitting procedure con- 
sists in transferring MC events from the CRs to 
the antiproton shadow and comparing the result 
with the data. 

To make such a comparison, we firstly adopted 
the following method. We obtained two kinds 
of Moon shadow, cast by all CRs and protons, 
respectively. After projecting them along the 
R.A. direction, we used a superposition of several 
Gaussian functions to describe the deficit event 
distribution in each shadow. Four Gaussian func- 
tions were found to be adequate for fitting both 
distributions within 5° from the Moon disc cen- 
ter. Let us name 9 the angular distance from the 
Moon disc center and f m (9) the Gaussian func- 
tion superposition describing the CR shadow. Let 
J~ P (0) be the proton shadow, obtained by impos- 
ing a given power law spectrum. The observed 
Moon shadow should be expressed by the follow- 
ing function: 

f moon (9) = (1 - r) f m {9) + rT ¥ {9) 
= (l-r)f m (9)+rT p (-9) 

(0 < r < 1) where the first term represents the 
deficit in CRs and the second term represents the 
deficit in antiprotons. This function must be fit- 
ted to the data to obtain the best value of r. 

We also applied a second method to determine 
the antiproton content in the cosmic radiation. 
Without introducing functions to parameterize 
the expectations, we directly compared the MC 
signal with the data. We performed a Maximum 
Likelihood fit using the p content as a free param- 
eter with the following procedure: 

I. the Moon shadow R.A. projection has been 



drawn both for data and MC. 

2. the MC Moon shadow has been split into 
a "matter" part plus an "antiproton" part, 
again so that the total amount of triggered 
events remains unchanged: 

§ MC {mat) — > $ MC (r;mat+p) = 
= (I - r)$ MC {mat) + r$ M c(p) 

3. for each antiproton to matter ratio, the 
expected Moon shadow R.A. projection 
®Mc( r i mat + P) is compared with the ex- 
perimental one via the calculation of the 
likelihood function: 

B 

logC{r) = Yl Niln[Ei(r)} - E^r) - ln(W) 
»=i 

where Nj is the number of experimental 
events included within the i-th bin, while 
Ei(r) is the number of events expected 
within the same bin, which is calculated by 
adding the contribution expected from MC 
(<&Mc{ r ', rll Q-t "Hp)) to the measured back- 
ground. 

Both methods described above give results con- 
sistent within 10%. The r parameter which best 
fits the expectations to the data turns out to be 
always negative, i.e. it assumes non-physical val- 
ues throughout the whole energy range investi- 
gated. With a direct comparison of the R.A. pro- 
jections, the r- values which maximize the like- 
lihood are: -0.076±0.040 and -0.144±0.085 for 
40<N<100 and N>100, respectively. The corre- 
sponding upper limits with 90% confidence level 
(c.L), according to the unified Fcldman & Cousins 
approach [16], are 0.034 and 0.041, respectively. 

Since the anti-shadow was assumed to be the 
mirror image of the proton shadow, we assume 
for the antiprotons the same median energy. The 
p/p ratio is $(p)/$(p) = l// p - (matter), 
therefore, being the assumed proton fraction 
/ p =73% for 40<N<100 and / p =71% for N>100 
[13], we obtain the following upper limits at 90% 
c.L: 0.05 for 40<N<100 and 0.06 for N>100. No- 
tice that the two values are similar, in spite of the 
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different multiplicity interval. It is a consequence 
of the combination of the two opposing effects of 
the angular resolution and of the geomagnetic de- 
viation. 

In Fig. H the ARGO-YBJ results are shown 
with all the available measurements. The solid 
curves refer to a theoretical calculations for a pure 
secondary production of antiprotons during the 
CR propagation in the Galaxy by Donato et al. 
[¥]. The curves was obtained using the appropri- 
ate solar modulation parameter for the PAMELA 
data taking period 2 . The long-dashed lines re- 
fer to a model of primary p production by anti- 
stars [5]. The rigidity-dependent confinement of 
CRs in the Galaxy is assumed oc R - * 5 , and the 
two curves represent the cases of S = 0.6, 0.7. 
The dot-dashed line refers to the contribution of 
p from the annihilation of a heavy dark matter 
particle [3]. The short-dashed line shows the cal- 
culation by Blasi and Serpico [T7] for secondary 
antiprotons including an additional p component 
produced and accelerated at CR sources. 




Figure 4. The p/p flux ratio obtained with the 
ARGO-YBJ experiment compared with all the 
available measurements and some theoretical cal- 
culations (see text). 



6. Conclusions 

The ARGO-YBJ experiment is observing the 
Moon shadow with high statistical significance at 
an energy threshold of a few hundred GeV. Using 
all data collected until November 2009, we set 
two upper limits on the p/p flux ratio: 5% at 
an energy of 1.4 TeV and 6% at 5 TeV with a 
confidence level of 90%. 

In the few-TeV range the ARGO-YBJ results 
are the lowest available, useful to constrain mod- 
els for antiproton production in antimatter do- 
mains. 
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